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ABSTRACT 



We report subarcsccond-rcsolution VLA imaging of four nearby spiral galax- 
ies: IC 342, Maffei II, NGC 2903, and NGC 6946. In each galaxy, compact radio 
continuum sources are identified in the central ~ 15" x 15" region. These com- 
pact sources are responsible for 20 - 30 % of the total emission from the central 
kpc of the host galaxies at 2 cm, but only ~5 - 10% at 6 cm. More than half 
of the compact sources appear to be Hll regions. The Hll regions with rising 
spectra must be fairly dense {rii ~ 10^ cm~^) and are presumably very young. 
The largest of these Hll regions require the excitation of 500 - 800 O stars, within 
regions of only few parsecs extent. These clusters approach the sizes expected for 
globular clusters. Thermal free-free emission from compact sources contributes 
more significantly at 2 cm, while diffuse synchrotron emission dominates at 6 
cm. The radio H ll regions are found near the centers of giant molecular clouds 
in projection, and do not have obvious visual counterparts. 

Subject headings: H II — regions galaxies: individual (IC 342, Maffei II, NGC 2903, 
NGC 6946) — galaxies: ISM — galaxies: starburst — galaxies: star clusters — 
galaxies: spiral — radio continuum: galaxies-infrared 
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1. Introduction 

Young massive star clusters, potential young globular clusters, have been revealed 
by ground-based imaging (Arp & Sandage 1985) and Hubble Space Telescope (HST) 
observations (see review by Whitmore 2003). These clusters appear to be a few Myr in age, 
10^ - lO'^ Mq in mass, and 10 pc or smaller in diameter. The youth of these clusters means 
that, unlike globular clusters, their luminosities are dominated by a population of massive 
stars. Young, massive star clusters appear to be common (Boker et al. 2002) in the centers 
of late-type spirals, and may play an important role in galactic nuclear activity. How these 
large star clusters form and the environments favoring their formation are not understood. 
The identification of young super star clusters in nearby galaxies gives us the chance to 
study these clusters in the process of formation. 

In nearby starburst galaxies, extremely compact sources with rising spectral indices 
at centimeter wavelengths have been discovered, which are thought to be "supernebulae" 
surrounding very young super star clusters (SSC) (Turner, Ho, & Beck 1998; Turner et 
al. 1999; Kobulnicky & Johnson 1999; Beck, Turner, & Kovo 2000; Tarchi et al. 2000; 
McDonald et al. 2002). These nebulae have the characteristics of dense, compact Hll 
regions in the Galaxy, and by analogy are also expected to be very young, ~1 Myr or less. 
However, they are much larger than Galactic compact H ll regions, requiring the excitation 
of hundreds to thousands of young O stars (Turner, Ho, & Beck 1998). These young 
clusters may be optically obscured by their natal gas and dust clouds (Wood & Churchwell 
1989), particularly in molecular gas-rich galactic centers. Free-free emission from embedded 
H II regions at radio wavelengths is unaffected by extinction from molecular clouds. Radio 
continuum emission can be used to search for the youngest, potentially embedded young 
H II regions in other galaxies. 

The goal of this investigation is to identify young SSC candidates in four nearby and 



-4- 



well-studied spiral galaxies, IC 342, Maffei II, NGC 2903, and NGC 6946 using subarcsecond 
radio continuum imaging. The centers of these galaxies are infrared-bright (Becklin et al. 
1980; Telesco & Harper 1980; Wynn-Wilhams & Becklin 1985; Ho et al. 1989) and molecular 
gas rich (Rickard et al. 1977; Young & Scoville 1982; Verter 1985). Previous, lower (~ 1") 
resolution VLA maps have shown that the four galaxies have concentrated regions of star 
formation in their centers (Turner & Ho 1983, 1994; Wynn- Williams & Becklin 1985). 
However, the low resolution of these maps, ~ 20 - 70 pc on the galaxy, combines emission 
from Hll regions, supernova remnants (SNRs), and extended non-thermal emission within 
the beam (Condon 1992). To isolate nebulae around individual star clusters we need 
subarcsecond resolution. 

Subarcsecond imaging with the VLA in its extended configuration maximizes 
sensitivity to bright and compact radio nebulae over low brightness disk synchrotron 
emission. Observing at shorter wavelengths, A < 2 cm, also minimizes the contribution of 
synchrotron emission, which falls with frequency. We present maps of the radio continuum 
emission at 6 and 2 cm in IC 342, Maffei II, NGC 2903, and NGC 6946 with the A 
configuration of the VLA. The images have resolutions of 0'.'3 at 6 cm and O'.'l at 2 cm, 
corresponding to size scales of 3 (IC 342 and Maffei 2) - 11 (NGC 6946) pc on the galaxies. 



2. Observations 

Table 1 lists basic properties for the four spiral galaxies. In all four there is prior 

evidence for compact radio continuum sources (van der Hulst, Crane, & Keel 1981, 1983; 
Turner & Ho 1983; Wynn-WiUiams & Beckhn 1985; Turner & Ho 1994). The radio data 
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were acquired at the NRAO Very Large Array^ using the A configuration. Previously 
unpubhshed 6 cm continuum data were obtained on 1988 November 18 and 20 (program 
AT98), and 2 cm data on 1983 October 12 (program ID AH141) and 1999 August 23 
(program AT227). The 2 cm observations in August 1999 were made using fast-switching, 
with cycle times of 2 minutes. On source integration times were ~ 1 hour at both 
wavelengths. Data calibration was done by using AIPS, following high frequency reduction 
procedures for the 2 cm cahbration. The sources 3C 48, 3C 138, 3C 147, and 3C 286 are 
used as fiux calibrators. The uncertainty in the absolute fiux scale is < 5 %. Observational 
parameters are shown in Table 2. 

To enhance the sensitivity, we combined our VLA A-configuration data with previously 
pubhshed VLA archival data hsted in Table 2. Data sets of A-, B-, or C-configurations 
with time on source > 10 minutes and phase center within 18 arcsec (1/10 of VLA primary 
beam at 2 cm) from centers of our measurements were used. Observations with B1950 
equinox coordinates were precessed to J2000 coordinates; uncertainty for this conversion is 
~20 mas. The uncertainty in the absolute position is determined by the B1950 calibrator 
positions, which are good to 50 mas. 

Spectral index measurements were done with a separate set of maps with matching 
{u,v) coverages at 6 and 2 cm, which include matched shortest baselines of A-configuration 
at 2 cm and longest basehnes of A-configuration at 6 cm. The largest angular scales sampled 
by the images are ~ 6". Fluxes and peak fiuxes are therefore lower hmits to the total fiux 
if extended emission is present. The images were then convolved to the same beamsize, to 
have matching maximum baseline lengths. Final rms noise levels in blank regions of maps 
are ~ 0.05, 0.04, 0.03, and 0.03 mJy/beam for IC 342, Maffei 2, NGC 2903, and NGC 6946 

^The National Radio Astronomy Observatory is a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated Universities, Inc. 



-6- 



respectively at 6 cm, and ~ 0.09, 0.1, 0.05, and 0.14 mJy/beam at 2 cm. 

While it is impossible to precisely match {u,v) coverages nor to quantify the 
differences without prior knowledge of the extended source structure, the {u,v) plane is 
well sampled and we estimate that differences between the maps caused by the differences 
in undersampled extended emission (> 10") must be less than the total resolved out flux 
divided by this area, which comes to less than our noise level. Since the compact sources 
we identify are all much less than 10" in extent, the effects on our fluxes with differences in 
{u,v) coverages will be undetectable. 

Interferometer maps are high-pass filtered images and it is difficult to assign 
uncertainties without knowing the extended structure, which is not known a priori. Point 
sources, and sources that are well-sampled by the (u,v) coverage will have uncertainties that 
are described well by the rms noise. Most of our sources fall into this category. We have 
fitted sizes to the extended sources to better estimate flux uncertainties for these sources. 

Although the differences between 6 and 2 cm spatial sampling may not affect compact 
source fluxes, we stress that there is extended unresolved flux in these maps. One generally 
compares single dish fluxes to mapped VLA fluxes to estimate how much extended emission 
is resolved out of the VLA maps. In this case, the disparity between single dish beams, 
typically a few arcminutes in size, and our sub-arcsecond beams makes this comparison 
difficult. We therefore compare our compact source fluxes to lower resolution VLA maps 
in B and C configuration (shown in Table 1) to see how much flux is resolved out; these 
lower resolution maps are sensitive to structures up to 45" across at 6 cm and 2 cm, except 
the case of NGC 2903, in which the C-array flux at 6 cm (sensitive to structure > 120" in 
size) has been quoted. The total mapped fluxes indicate that our maps recover 15-45% of 
the 6 cm flux in the lower resolution observations, and 20-45% of the 2 cm flux, in both 
compact and extended components. 
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3. The Radio Continuum Images and the Compact Radio Sources 

The 6 cm radio continuum maps and K-hand (2.2//m) images of the four galaxies are 
shown in Figure 1. The naturally- weighted radio contours are overlaid on 2 Micron All Sky 
Survey {2MASS) i^-band images on the left. In each of the galaxies the nucleus is a strong 
2 /im source and the 6 cm continuum is strongest at the center. We see emission from 
the central ~ 150 pc except in the case of Maffei II, in which the emission extends over a 
distance of ~ 350 pc. 

The 6 cm maps on the right-hand side of Figure 1, contoured at the same 4 a levels as 
in the overlay on the left-hand side, are closeups of the radio emission region within the 
central 16" - 26". In each source there are strong compact sources embedded in extended 
emission. The western part of the extended emission in IC 342 contains 5 sources aligned 
from north-west to south-east (named A - E in Table 3) while the eastern part has fewer 
sources and their emission is also significantly weaker. Ten sources are identified in Maffei 
II, spread out over 10" (~ 240 pc) north-south on the eastern edge of the extended emission. 
6 cm continuum emission in NGC 2903 and NGC 6949 is confined to single central sources 
surrounded by extended emission. The diffuse emission is "patchy" because the large scale 
extended emission has been resolved out by the high resolution. We discuss features in the 
individual galaxies in subsequent sections. Peak flux densities at 6 cm are 1.2, 2.6, 0.4, and 
1.9 mJy/beam for IC 342, Maffei II, NGC 2903, and NGC 6946, respectively. 

The radio maps are overlaid on HST F-band and Ha images in Figures 2-5. For 
convenience in comparing to the 2 cm, the 6 cm maps in Figures 2-5 were restricted to 
the same {u,v) coverage (~ 30 - 600 kX), and convolved to the same beamsize. The radio 
and the HST images have similar resolutions. The uncertainty in the registration of the 
images is ~1", hmited by the astrometry of the HST images. In IC 342, Maffei 2, and NGC 
2903 the radio sources tend to lie in regions of apparent dust obscuration. In NGC 6946 the 
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extended central radio continuum source nearly aligns with a bright optical source at the 
center. The uncertainty in registration does not allow us to determine how precisely they 
align. 

In Table 3, we list compact sources detected in the central 15" radius of each galaxy 
and their fluxes at 6 cm and 2 cm, with proper {u,v) restrictions described in §2. The 
sources meet at least one of the following criteria: (1) 5 cr detection at the peak intensity 
at one wavelength, or (2) 4 a emission detection in both wavebands. The sources are 
named alphabetically following radio convention in order of their peak flux densities at 
6 cm. Positions and deconvolved sizes of sources were obtained from the 2 cm data fit 
assuming a single Gaussian source using the IMFIT task in AIPS. For sources absent at 2 
cm, Gaussians were fit to the 6 cm maps. The uncertainty in size is < 30 % of the beam. 
The total integrated fiux densities were measured using the AIPS task IRING. The upper 
hmits quoted for non-detections are 3 a. The systematic error of fiux is expected to be < 
5%, related to the uncertainties in the absolute fiux cahbration fiuxes. We separated this 
uncertainty from signal to noise, which is listed in the Table 3. 

In Table 4, we list the properties of each compact radio source. The spectral index 
is defined as oc i/". The uncertainties of a6-2 5 the spectral index between 2 and 
6 cm wavelengths, are derived with -\/2 a variances of fiux densities in both wavebands; 
uncertainties are typically ±0.2-0.3. For sources A - D and I in the extended emission 
complex of NGC 6946, we derive the a from their intensities, which are well-determined, 
rather than their fiuxes, which are more uncertain. Based on the values of q;6-2 (discussed in 
the next section), we classified sources as optically thick Hll regions (Hll, thick), optically 
thin Hll regions (HII, thin), supernova remnants (SNR), or radio supernovae (RSN). Source 
types are listed in the fourth and fifth columns in Table 4. For the cases identified as 
H II regions, A^l^/c in the sixth column is the Lyman continuum rate derived from the fiux 
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density at 2 cm on the assumption of optically thin emission and electron temperature in 
H II region, ~ 10,000 K. If the emission is optically thick, this value will be a lower 
limit to the true A^l^c- Electron temperatures of H ll regions tend to be cooler toward the 
high-metallicity centers of spiral galaxies (Pagel et al. 1979; Sedwick & AUer 1981). The 
values of A^Lyc will be ~ 15 % lower than listed, if Tg — 6,000 K rather than the adopted 
10"^ K. No7 is the number of standard 07 stars needed to generate N^yc- We also assume 
ionization-bonded nebulae; if UV photons escape the H ll region, our OB luminosities will 
be underestimates of the true Lyman continuum rates. 

4. Discussion 

4.1. Identification of Luminous Compact Hii Regions from their Radio 

Spectra 

The radio images reveal a population of bright compact radio sources in the centers of 
IC 342, Maffei 2, NGC 6946, and NGC 2903. The sources are comparable in extent to the 
0'.'3 beam, or ~ 5 - 10 pc in diameter. What are these sources? 

The cm-wave continuum emission from nearby normal galaxies is a mixture of 
synchrotron emission from SNRs and disk cosmic rays, and thermal bremsstrahlung 
emission from H ll regions. Spectral index a is an indicator of the emission mechanism: 
optically thin free- free emission has a ~ —0.1, while optically thin non-thermal emission by 
synchrotron radiation will have steeper negative spectral indices. Galactic SNR have ant ~ 
-0.6 to -0.2, with the majority in the range -0.5 to -0.4 (Green 1984; Mills, Turtle, Little, 
& Durdin 1984; Weiler et al. 1986). Steeper 0:6-2 < —1.0 are characteristic of short-lived 
RSNe (Weiler et al. 1986), and of extended synchrotron disk emission. 

Hll regions are optically thick at wavelengths longer than a turnover frequency 



-10- 



that depends on the electron temperature and emission measure, EM= J dl. At long 
wavelengths, Hll regions are bright, and their fluxes rise with frequency until the turnover 
frequency. The higher the emission measure, the higher the turnover frequency. Spectra of 
diffuse or "classical" (EM < 10^ cm~^ pc) Hll regions that appear in visible wavelengths 
typically turn over at meter wavelengths and are optically thin at cm wavelengths. A 
cm- wave turnover frequency is characteristic of compact or ultracompact Hll regions 
(Mezger & Henderson 1967), which have EM > 10® — 10^ cm~^ pc. 

Beyond the Local Group, ultracompact Hll regions typical of those in our Galactic 
disk, excited by individual O stars or small clusters, are small in size (Habing & Israel 
1979) and thus too faint to be seen (Turner & Ho 1994). The compact component (Mezger 
et al. 1967) of the luminous Galactic Hll region W49 would fail to meet our detection 
criteria in our closest galaxy (IC 342), as would the radio source (de Pree et al. 1999) 
of NGC 3606. However, if the clusters are somewhat larger, with at least 200 O stars 
surrounded by an ionization-bounded volume of gas, even dense, compact H ll regions, with 
EM > 10® — 10^ cm~^ pc, can be detected in external galaxies (Turner, Ho, & Beck 1998). 
Even for the largest clusters (A^oy = 10"^), dense, compact Hll regions will still be small, 
less than 2 - 3 pc in size. Since large numbers of massive stars can also be accompanied by 
supernovae, thermal free-free emission from these H ll regions can be mixed with nonthermal 
emission within a single beam if the resolution is low (Turner, Ho, & Beck 1998; Turner, 
Beck, & Ho 2000). Confusion is minimized if the resolution is high, and the galaxy is 
nearby, so that the beam closely matches the size of a single H ll region. Subarcsecond 
imaging maximizes our sensitivity to these sources. 

More than half of the 38 identifled compact sources of Table 4 appear to be H ll 
regions. A significant fraction may have rising spectra, indicating that they are at least 
partly optically thick at 6 cm. Given the likely presence of nonthermal synchrotron emission 
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in the vicinity of these star-forming regions, the spectral indices could systematically be 
even higher than we measure here, due to spectral confusion (Turner, Beck, & Ho 2000). 
Rising cm-wave spectra imply that the rms electron densities are ~ 10^ — 10"^ cm~^. That a 
large percentage of our sources are ultracompact H ll regions might be expected from the 
subarcsecond resolution, which picks out these high brightness regions. There may be lower 
density, extended and less bright H ll regions that are below the brightness limit of these 
maps. 



4.2. Compact Nebulae as Young Star Formation Tracers 

4-2.1. Properties of the Compact Hll Regions 

The sizes of the compact H ll regions listed in Table 3 are obtained from deconvolution 
of the beam from the observed source in the 2 cm map, assuming that the sources are 
gaussian. (They are probably not gaussian in structure, but this gives us approximate source 
sizes.) Uncertainties in the sizes, due to beam size, source structure, and signal-to-noise 
ratio, are < 30 % of the beam. The derived source sizes in IC 342, Maffei II, and NGC 2903 
are ~ 0.3", but ~0.5" in NGC 6946. Sources with larger deconvolved sizes in NGC 6946 are 
nebulae joined with the large emission complex at the center, and suffer from confusion. 
Seeing may be a problem at 2 cm. The similarity of the source sizes suggest this may be the 
case; it is likely, therefore, that the sizes we obtain are upper limits to the true source sizes. 

The deconvolved sizes correspond to diameters of ~ 6 - 14 pc. These H ll regions are 
significantly smaller than the 30 Doradus H ll region (~ 200 pc in diameter for region of 
EM > 10^ pc cm~^; see Mills et al. 1978 and Kennicutt 1984) in the Large Magellanic 
Cloud (LMC). They are larger but similar in nature to the Galactic compact Hll regions in 
W49 (Mezger et al. 1967; Conti & Blum 2002) and NGC 3603 (de Free et al. 1999; Sung 
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& Bessell 2004; Miicke et al. 2002). One possibility for the fact that these Hll regions 
are relatively compact and dense is that they are younger than 30 Doradus or NGC 3603. 
Because these regions are nearly unresolved, and because we do not have good radio SEDs, 
at this stage we cannot model their properties with precision. 

The estimated sizes of these nebulae are similar to the pc-sized diameters of the 
supernebulae found in NGC 5253 and II Zw 40 (Turner & Beck 2004; Beck et al. 2002; 
Meier et al. 2002). These Hll regions have ionized gas densities at least as high as 
Tii ~ 10^ — 10^ cm~^; they are dense, with high pressures. The compactness of these 
optically thick H ll regions implies some degree of youth, although H ll regions around 
massive young clusters may differ dynamically and evolve differently from Galactic compact 
Hll regions (Kroupa & Boily 2002; Turner et al. 2003; Turner & Beck 2004). 

The radio fluxes of these compact H ll regions require that they harbor hundreds of 
massive stars. The luminosities of the clusters exciting the H ll regions can be inferred 
from their Lyman continuum rates using Lqb = 2 x 10~^^ ^Lyc for a Salpeter IMF. The 
inferred luminosities are not very sensitive to IMF since luminosity, like ionizing photons, is 
dominated by the most massive stars. The cluster luminosities inferred from the free-free 
emission are ~ 10^ — 10^ Lq. Based on their luminosities, the largest of these radio clusters 
have roughly an order of magnitude more stars than the R 136 cluster within 30 Doradus 
(Massey & Hunter 1998; Zinnecker et al. 2002), or in the clusters within the large Galactic 
Hll region NGC 3603 (Eisenhauer et al. 1998; de Free et al. 1999) or W49 (Conti & Blum 
2002). 

The total Lyman continuum rate, N^yc, is derived from the total flux from thermal 
compact sources, assuming the emission is optically thin and the Hll regions are 
ionization-bounded. These values are compared with NLyc rates derived from millimeter 
continuum, which are unlikley to suffer from opacity effects. In IC 342, the estimated A^l^c 
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is ~ 6 X 10^^ for six thermal compact sources. This can be compared with the value 
of Niyc ~ 3 X 10^^ obtained by Meier & Turner (2001; corrected for distance) for the 
central 16" region derived from 2.6 mm continuum, and N^yc ~ 1.2 x 10^^ for the two 
major 2.6 cm peaks corresponding to the central starburst regions. Over the central 7" x 2" 
region of NGC 6946, NLyc is 3 ± 0.2 x 10^^ as derived from the 2.7 mm images of Meier 
& Turner (2004; corrected to our distance). Similar agreement is also found in Maffei II 
between the A^l^c inferred by thermal fluxes in mm (Meier, Turner, & Hurt 2006) and cm 
(this work). These numbers imply that the fluxes of compact sources make up at least 50 % 
of the total Niyc photons seen towards the starburst peaks in lower resolution images. The 
fact that the mm emission gives higher values for Ni^yc suggests that some of the thermal 
sources are optically thick at cm wavelengths, or that some of the thermal flux is extended, 
or both. 

4-2.2. Compact Nebulae and Molecular Clouds; Possible Star Formation Triggers 

The existence of massive stars within these clusters indicates their youth. Because the 
local star formation is expected to be suppressed by massive star feedback of mechanical 
energy through winds and supernovae, we expect the stars in these large clusters to be 
coeval to within a few Myr. Thus they are presumably located inside and co-moving with 
their parent molecular clouds, and should be close to the sites at which they formed. 
Given that these Hll regions are compact and quite dense, it is likely that they may be 
significantly younger than the ~ 10 Myr we might infer from the presence of massive stars. 

The young, massive clusters in IC 342 and Maffei II (shown in Figure 2 and 3) are 
aligned with dust lanes and molecular gas. These nebulae are found in coherent structures 
extending over regions 110 pc (IC 342) - 240 pc (Maffei 2) in size. The sound crossing 
times for these distances are 10 - 25 Myr, if the sound speed is ~10 km/s. This timescale 
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is comparable to or longer than the lifetimes of massive stars. Prom the youth inferred 
from the high gas densities and from the spatial distribution of these H ll regions along the 
CO arms, we conclude that the formation of young, massive stellar clusters in compact 
nebulae found in IC 342 and Maffei II is mainly triggered by large scale resonant features 
(i.e., gas arms associated with spiral density wave/bar features; Lo et al. 1984; Ishizuki et 
al. 1990; Ishiguro et al. 1989; Hurt & Turner 1991; Schinnerer, Boker, & Meier 2003; Meier 
& Turner 2004, 2005). This is consistent with their locations with respect to optical dust 
lanes (Figures 2-5). 

Comparing with the HST images (Figure 2 - 5), we find that only a small fraction of 
the radio H ll regions have possible optical counterparts. The extinction in embedded H ll 
regions is often internal to the nebulae (Beck et al. 1986; Kawara et al. 1989; Ho, Beck, 
& Turner 1990), but it is also possible that the compact nebulae are hidden behind the 
foreground clouds or embedded in their host giant molecular clouds. If we assume that all 
of the luminosity comes out in the infrared, the implied Ljr of our compact sources amounts 
to only 3 - 5 % of the total IR luminosities of the galaxies (Table 1). 

The young nebulae are all located within the central ~ 100 pc except in Maffei II, 
where the H II regions extend over ~ 250 pc. To support the formation of numerous massive 
stars, there should be high molecular gas concentrations (Young & Scoville 1991). In all 
four galaxies, the central regions have molecular gas surface densities 3-9 times higher 
(Hurt & Turner 1991; Turner & Hurt 1992; Meier & Turner 2001, 2004) than do typical 
spiral galaxies on size scales of 6" (Heifer et al. 2003). The high concentrations of molecular 
gas {Ugas = 600 - 2000 pc''^ within central 6"; Hurt & Turner 1991; Heifer et al. 2003) 
lead to high internal visual extinctions Ay of up to tens to hundreds of magnitudes (Ho, 
Beck, & Turner 1990; Hurt & Turner 1991; Turner & Hurt 1992) toward these galactic 
centers. 
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The young star clusters associated with the radio nebulae are preferentially found 
projected onto the centers of giant molecular clouds (§4.3). This means that there is a 
wealth of gas around the star clusters to fuel more star formation, unlike the high star 
formation efficiency inferred from the lack of gas in the super star clusters found in the 
dwarf galaxy NGC 5253 (Turner et al. 1997; Meier et al. 2002; Turner et al. 2002). Given 
the proximity of molecular clouds to these clusters, and their location in the centers of 
large spiral galaxies, it is unlikely that any of the young clusters associated with the radio 
nebulae in these four spirals will become long-lived, gravitationally bound clusters, i.e., 
globular clusters. 

4.3. Upper Part of the Luminosity P\inction of Young H ll Regions 

Nearly two dozen radio sources in the four galaxies were found to have flat or rising 
spectra (Table 4). These are presumably thermal nebulae excited by massive stars, for 
which the 2 cm fluxes give a lower limit to the required Niy,, (§ 4.1). In Figure 6, we show 
the Lyman continuum luminosity function (LF) of these H ll regions. The corresponding 
Ha luminosities have been calculated from the A^Lyc assuming no extinction, for comparison 
with other studies. The sample is small, and the statistical error bars are relatively large. 
However there are some trends that appear to be real and should be studied with a larger 
sample. 

We have fit the LF with a power law, and also with a two component power law. 
The two component fit is significantly better than the single component fit. The LF has 
a steeper slope for nebulae with Ha luminosity > 2 x 10^^ erg sec~^ (corresponding to 
Lob = 10^'^ Lq, about 150 07 stars). The solid fine in Figure 6 represents power law 
models {N{L) oc L"-dL) of LF with index of —2.2. The power law index we obtain is close 
to the average value of —2.0 in 30 nearby spiral and irregular galaxies from an optical 
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study of Hll regions by Kennicutt, Edgar, & Hodge (1989). However, the slope of the LF 
for the more luminous Hll regions of our sample is shallower than the slope in the two 
component power law LF, or "type 11" LF, found in optical studies of other spiral galaxies, 
M 51, NGC 3521, NGC 3627, NGC 4736 (Kennicutt, Edgar, & Hodge 1989), and NGC 
7714 ( Gonzalez- Delgado et al. 1995). The optical studies also show the turnover point at 
lower luminosity (at L^^ ~ 10^®'^ — 10^^). A possible explanation for this difference might 
be that the luminous sources we found are very young and obscured, thus not detected in 
the optical, although this may also be small-number statistics. 

The flattening of the LF at the faint end may be caused by incompleteness due to 
both hmited sensitivity and partially resolving out extended H ll regions. We cannot detect 
unresolved (diameter < 15 pc) nebulae with NLyc < 4 x 10™ (3cr), corresponding to 
40 O stars in IC 342 and 120 stars in NGC 6946. However, we also cannot detect more 
luminous, but lower brightness regions if they are large and resolved (size ^ 15 pc); in this 
limit we are sensitive to Hll regions with EM > 10^ cm~^ pc. 

We did the same thermal source identification and luminosity calculation for radio 
sources in M 82, another nearby star burst galaxy, based on radio studies of Allen (1999). 
Allen's source list is more complete, with Ni^y^ < 2 x 10"'''^ s^^ (20 07 stars) and radio 
nebulae with size < 50 pc. The LF of the 21 sources in M82 with flat or rising spectra 
between 2 and 6 cm shows very similar features in slope and the break point in LF to what 
we find in the four spirals (Tsai et al. 2006) . The agreement of LF shape of M 82 and 
our sample suggests that the broken power law may be real. Larger samples would help to 
clarify this issue. 

The steep power law LF suggests a cutoff at the high luminosity end in the young 
clusters associated with radio nebulae in spiral galaxies. The upper end of the luminosity 
function is around Lh„ ~ 10^*^ erg s~^. This corresponds to ~ 750 07 stars. We observe 
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two such clusters, in Maffei 2 and NGC 2903, and their nebulae are < 5 pc in extent. The 
total mass of such a cluster would be ~ 2.2 xlO^ Mq (assuming Salpeter IMF), consistent 
with the masses of Galactic globular clusters. 

5. Individual Galaxies 
5.1. IC 342 

IC 342 is a face-on giant Scd spiral (Abies 1971; Crosthwaite et al. 2000, 2001) close to 
the Galactic plane. Its distance is uncertain: Saha, Claver, & Hoessel (2002) used the P - 
L relation for Cepheids and suggest 3.3 Mpc, close to the 3.9 Mpc given by Tully & Fisher 
(1988), but earlier measurements found ~ 2 Mpc in optical spectroscopic and photometric 
studies (McCall 1989; Karachentsev & Tikhonov 1993). Here we adopt 3.3 Mpc (1" = 16 
pc). 

The core of IC 342 contains a central molecular ring 5" - 6" in radius and a nuclear CO 
disk ~ 3.5" (~ 55 pc) in diameter (Schinnerer, Boker, & Meier 2003). The 5" - 6" nuclear 
disk contains ~ 2 x 10^ Mq molecular gas and the central 100 pc contains ~ 5 x 10^ Mq of 
molecular gas (Meier & Turner 2001). The nuclear star cluster (NSC) of IC 342 has been 
studied in optical and infrared wavelengths using HST and IRTF. The NSC of IC 342 has 
mass M Ri 6 X 10^ Mq and age IQ^-^"'''-^ years. The upper limit on a central black hole is 
5 X 10^ Mq (Boker, van der Marel, & Vacca 1999). Outside the central star cluster, there is 
current star formation indicated by bright IR and radio emission clumps, each containing 
10^ - 10^ OB stars, and the total number of equivalent O stars is about 5,000 (Becklin et al. 
1980; Turner & Ho 1983). Our measurements resolve out most of the extended component 
of emission so our calculated total O star population, which detects compact sources only, 
is over an order of magnitude less. 
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The radio continuum images (Figure 2) reveal twelve compact radio continuum sources 
on the two spiral arms and the central molecular ring of IC 342. 10 % of the total radio 
emission at 6 cm and 20 % at 2 cm for the central kpc arises in compact sources. Inferring 
the infrared luminosities from their radio fluxes and N^yc {"Lob"), we find that the compact 
H II regions in the nucleus contribute only ~ 5 % of the total infrared luminosity of IC 342. 
2 cm sources A, C, E, G, and I - K are distributed along the inner edge of the dust lane and 
the western molecular spiral arm near the nucleus (Schinnerer, Boker, & Meier 2003). The 
western arm of the central molecular ring clearly has more energetic star formation than the 
eastern arm, as is evident from lower resolution observations. The star formation in IC 342 
appears to be spatially distributed in the nuclear region, although somewhat concentrated 
to the west near the northern end of the western molecular arm (Meier & Turner 2005). 

We identify a number of likely SNR. Source G in the eastern central molecular arm has 
a ~ —0.3. Four sources, B, D, F, and H, are detected only at 6 cm. Sources A and B have 
a ~ —0.5 and —0.6 respectively, close to the conventional a of Galactic SNR. The steep 
spectrum of Source D makes it a candidate for a RSN, which typically decline rapidly over 
the span of a few years. However, no obvious change in the 6 cm flux of source D is found 
in the data of IC 342 we have collected over multiple epochs separated over 15 years. Its 
diameter (~ 0.6" — 0.4", or 10 pc to 6 pc) and luminosity Lg^n = 2.9 x lO^^erg sec"^ Hz~^ 
can be explained better by a SNR. The unusual negative spectral index of source D might 
be due to the flux contamination by nearby non-thermal source F. Sources F and H show 
no strong emission at 2 cm, and have upper limits to their spectral indices of a < —0.5. 
None of these sources. A, B, D, or F - H, arc associated with obvious visible star clusters. 
As SNR candidates, they are the likely end stages of massive stars near the young star 
cluster systems. 

We identify flve compact sources that appear to be compact Hll regions. Sources 
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C, I, J, and K in the north-west and source L in the east have spectra consistent with 
optically thick free-free emission. If the emission were optically thin, the excitation of these 
compact nebulae requires clusters of at least 70 - 210 07 stars each. These are lower limits 
to the true cluster population, because the emission is optically thick, and because of the 
possibihty of dust absorption of the UV photons. The inferred cluster masses are then 
2.1 - 6.3 X 10^ Mq. These are large clusters comparable to or a few times larger than 30 
Doradus in the LMC, but much more compact. These nebulae are less than 8 pc in size as 
compared to 200 pc for 30 Doradus, and hence probably younger. These star clusters are 
not detected in the HST image, implying that Ay > 5 mag, or Nh ^ 1.0 x 10^^ cm~^, if 
the ZAMS star clusters with Salpeter IMF are assumed. The 4" resolution 10 /xm image of 
IC 342 (Beckhn et al. 1980) shows two peaks at (~ 03'*46"*48.3^; +68°05'47") and (~ 47.6^; 
46"), close to the radio Hll regions we identify here. The visual extinction estimated from 
silicate absorption feature is ~ 15 mag toward these two 10 /xm peaks (Becklin et al. 1980), 
consistent with our argument of high extinction toward these H ll regions. 

To the west. Source E has a flatter radio spectrum and appears to be an optically thin 
source. This Hll region requires about 70 O stars, comparable to 30 Doradus. A ZAMS 
star cluster with 70 O stars in IC 342 should have apparent brightness my = 13.8 when 
Galactic extinction of Ay = 1.85 {Nh ~ 3.5 x 10^^ cm~^) toward IC 342 is applied; these 
clusters do not appear in the HST optical image (Fig 3, also see Boker, van der Marel, & 
Vacca 1999). 

We do not detect significant radio continuum emission at small size scales from the 
nuclear star cluster to a 3 cr hmit of 6.5 x 10^°sec~^, or 70 equivalent 07 stars. This is 
consistent with a somewhat older stellar population there (60 Myr Boker, van der Marel, 
& Vacca 1999). No significant active galactic nucleus has been detected with these high 
resolution measurements, but our own Galactic center at a similar resolution (Brown & 
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Johnston 1983) would be below our detection limit at this distance (~ 0.1 mJy at 2 cm). 

5.2. Maffei II 

Maffei II, discovered in the infrared by Maffei (1968), is a spiral galaxy with a strong 
and long bar evident in near-infrared images (Hurt, Merrill, Gatley, & Turner 1993). It 
is located behind the Galactic plane {I — 136.50°, b — —0.33°) at Ay ~ 5 - 10 mag of 
Galactic extinction. The nucleus of the galaxy is actively forming stars (Rickard & Harvey 
1983; Ho, Beck, & Turner 1990; Telesco, Dressel, & Wolstencroft 1993; Seaquist, Pfund, & 
Bignell 1976; Turner & Ho 1994). The distance of Maffei II was recently set at 2.8 Mpc by 
Karachentsev, Sharina, Dolphin, & Grebel (2003) using the TuUy-Fisher relation; however 
the HI linewidth suffers from Galactic absorption, particularly in lower resolution spectra 
(Hurt, Turner, & Ho 1996) so the actual Hi linewidth is larger than what they adopt. In 
this paper, we use the distance of 5 Mpc estimated from the diameters of the largest H ll 
regions in Maffei II (Spinrad et al. 1973) and by the TuUy-Fisher analysis of Hurt, Turner, 
Ho, & Martin (1993). The uncertainty in the distance to Maffei 2 introduces a factor of 
3 uncertainty in estimates of massive star numbers and total masses, and a factor of 2 
uncertainty in source sizes. 

In Figure 3 we show the radio continuum maps of Maffei II. Ten sources are identified 
in both 2 and 6 cm maps and sources B and J are identified only at 6 cm (see Table 3 and 
Fig 3). The compact sources he on a twisted, inverse "S" shaped line north-south, following 
the trend of larger scale gas distribution and coinciding with the dust lane shown in the 
HST F-band image. The radio continuum emission is located along the elongated bar of 
45" length seen in ^^CO (Ishiguro et al. 1989) and ^^CO (Hurt & Turner 1991) images. 
Northern sources B, C, and F are associated with giant molecular cloud (GMC) Dl of 
Meier, Turner, & Hurt (2006). Sources A, D, H, I, and J are close to the center of GMC 
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E. The southern sources, E and G, are between GMC E and F. Strong detections in both 6 
and 2 cm provide good spectral indices and sizes for the radio sources in Maffei II, except 
for source B and J which have no significant detection at 2 cm. 

The total flux densities from compact sources constitute ~15 % of the total emission 
of the central kpc of Maffei II at 6 cm and ~30% at 2 cm (for structures < 45"; see §2). 
Turner & Ho (1994), using the VLA in B array at 6 cm and C array at 2 cm, found a total 
number of O stars in good agreement with our results. The star formation in the center of 
Maffei 2 is concentrated in large, compact clusters aligned in a linear structure along the 
northern CO arm. The predicted infrared luminosity Ljji — Lqb ~ 5.3 x 10^ Lq from these 
point sources is ~5 % of the total Lm of Maffei II. Except for three SNR/RSN candidates, 
the compact radio continuum sources we detect are thermal free- free emission from H ll 
regions. Most of the thermal sources in Maffei II have relatively flat spectra, compared with 
rising spectrum sources in NGC 2903 and NGC 6946. 

Spectral indices a ~ —0.1 for sources A, D and I suggest that these are star clusters 
with optically thin thermal free-free radiation. All three are in the center of the radio 
emission near GMC Dl. These Hll regions of diameter < 10 pc each require ~ 200 - 600 
exciting O stars. The inferred masses of the clusters associated with sources C and G are 
~6-18xl0^ Mq, slightly larger than the young massive star clusters in our Galactic 
Center. 

Sources C, E, G, and H have rising spectra indicating that they are compact Hll 
regions. These Hll regions are brighter and more luminous than the ones in IC 342. The 
lower limit of A^lj/c = 2.9 - 4.9 x 10^^ requires clusters of more than 300 - 500 07 
stars. The inferred cluster masses are Mq ~ 1 - 1.5 x 10^ Mq. The mass of source C, 
Mc ~ 1.5 X 10^ Mq is only a smaU fraction of the 26 x 10^ Mq upper hmit mass of GMC C 
(Meier, Turner, & Hurt 2006). The star formation efficiency therefore appears to be of order 
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one percent, which is similar to the Galactic value on these scales. Prom the C^^0(2-l) 
map of Meier, Turner, & Hurt (2006), Source E and G seem to be at the gas bridge in the 
middle of GMC Dl and GMC E3. The predicted visual brightness of these rising spectrum 
sources with the assumption of a Salpeter IMF is my ~ 25 mag with Galactic extinction, 
so these sources are too faint and highly obscured to be detected via HST observations. 

Source F is probably a SNR, with a — 0.4. The luminosity of source F at 6 cm is 
L&cm ~ 4.6 X 10^^ erg sec~^ Hz~^ and the deconvolved size is 11 x 8 pc, characteristic 
of a young Type II SNR. Sources B and J show no significant counterpart at 2 cm, and 
their negative spectral indexes suggest that they are RSN and SNR, respectively. All three 
non-thermal sources, F, B, and J, are within > 1" (25 pc in projection) to the nearby 
optically thick thermal sources, C and H. The locations of these SNR/RSN candidates 
suggest they might be associated with the compact H ll regions. 

5.3. NGC 2903 

NGC 2903 is a large and bright "hot spot" spiral, with an unusually complex structure 
of both early- and late-type stars and Hll regions at its center (Oka et al. 1974; Sersic & 
Pastoriza 1967; Morgan 1958; Sersic 1973). It is at a distance of 8.9 Mpc (Drozdovsky & 
Karachentsev 2000). Wynn- Williams & Becklin (1985) studied the "hot-spot" galaxy in 
the radio and infrared, and found that its radio emission is predominantly nonthermal. 
However, in some regions the free-free emission contributes substantially at 2 cm. NGC 
2903 has a well-studied star-forming ring. IR maps find Hll regions in the ring with 
luminosities approaching that of 30 Doradus (Alonso-Herrero, Ryder, & Knapen 2001). IR 
spectroscopy gives an age range of 4 - 7 xlO^ yr for the embedded star clusters while the 
Hll regions are thought to be only 1-4 xlO^ yr of age (Alonso-Herrero, Ryder, & Knapen 
2001). Molecular gas in NGC 2903 is concentrated in the central few tens of arcseconds 
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(Sheth et al. 2002; Heifer et al. 2003). 

In Figure 4 we show NGC 2903 in radio continuum maps and HST WFPC2 images 
( F-band and Ha). About 20% of the total flux density at 2 cm emission in the central kpc 
comes from compact sources B, C, E, F, and G. At 6 cm ~5 % is due to sources A - E. This 
agrees with Wynn- Williams & Becklin (1985) on the contribution of compact sources at 
radio wavelengths. The northern peak near 09'*32"*09.99*, 21°30'12.7" seen in low resolution 
2 cm maps (Wynn- Williams & Beckhn 1985) does not have a corresponding compact source 
in our images, and therefore must be extended over >2" scales. 

K and if-band NICMOS and Pa a images of the nucleus of NGC 2903 (Alonso-Herrero, 
Ryder, & Knapen 2001) show an absence of ionized gas, consistent with what is seen in 
our 2 cm and 6 cm maps. The nucleus must consist of a relatively old stellar population 
dominated by giants and supergiant stars. The northern peak of the lower resolution 2 cm 
maps of Wynn- Williams & Beckhn (1985) corresponds to compact Source B. It has a rising 
spectrum of o; ~ 0.8 through cm wavelengths. Excitation of Source B requires at least 
800 07 stars if the free-free emission is thin and more if thick. This result is ~ 2.6 times 
higher than N^yc derived by Alonso-Herrero, Ryder, & Knapen (2001). It suggests visual 
extinctions of 2.7 mag higher than what Alonso-Herrero, Ryder, & Knapen (2001) obtained, 
thus Ay = 6.5, assumed if ^v/Ai.g/^m apphed (CardeUi, Clayton, & Mathis 1989). This 
difference could be explained if there is significant extinction at Pa a. The brightness 
temperature of Source B indicates a size of 0.025", or 1.1 pc for a nebula with T^. =10,000 
K. Source B is a good candidate for a super star cluster nebula. 

Sources A, D, E, and G are within the extended emission east of the northern peak 
(Wynn- Williams & Becklin 1985) near Source B. All of them have counterparts in Paa 
(Alonso-Herrero, Ryder, & Knapen 2001) except Source E. Source E is close to source "b" 
of Oka et al. (1974), which they claim to be a free-free emitter. Source E also appears 
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coincident with cluster 2 in the HST near-infrared images Alonso-Herrero, Ryder, & Knapen 
(2001) and the AO if '-band image (Perez- Ramirez et al. 2000). Sources A, B, D, and G 
correspond to sources HI, H4, H2, and H3 of Perez-Ramirez et al. (2000), respectively. 
In our maps. Source G has been identified only at 2 cm, and shows a very steep rising 
spectrum. Its spectral index, a > 1.0, implies that this is a possible compact Hll region 
with a high electron density and should therefore be very young. Source D, detected at 
6 cm but not at 2 cm,is probably a SNR. Source A has a flat slope of o; ~ O.Otoi. It could 
be a small Hll region, with ~ 210 equivalent O stars, or, within the uncertainties, it could 
also be a single young Crab-hke SNR with the luminosity of Gas A. 

Southern Sources G and F have thermal spectra. Source G coincides with the southern 
peak in the low resolution 2 cm map of Wynn-WiUiams & Beckhn (1985). The flux density 
at 2 cm is about a quarter the peak in the G-configuration maps, so this source is extended. 
Source G contains about 230 O stars in 0.4" x 0.3", based on its 2 cm flux, and 100 
equivalent O stars based on the Pa a measurement of Alonso-Herrero, Ryder, & Knapen 
(2001); the discrepancy is probably due to extinction at Pa«. Source F is an optically thick 
Hll region with spectral index a > 1.2. The cluster inferred to be exciting source F has 
Lob ~ 1-6 x lO^L© and mass Mc < 2.3 x lO^M© which is about the mass of a globular 
cluster. 



5.4. NGC 6946 

NGC 6946 is a barred spiral galaxy close to the Galactic plane (/ = 95.72°, b = 11.67°). 
It is a bright mid-infrared, (Rieke & Lebofsky 1978), GO (Ball et al. 1985; Sofue et al. 1988; 
Wehachew, Gasoli, & Gombes 1988; Gasoh, Dupraz, & Gombes 1992; Regan & Vogel 1995; 
Sakamoto, Okumura, Ishizuki, & Scoville 1999), and radio (Turner & Ho 1983) source. 
NGG 6946 contains a moderate starburst of Ljji ~ 2.2 x lO^L© and age range 7-20 Myr 



-25- 



(Engelbracht, Rieke, Rieke, & Latter 1996), and 1 x 10^ M© of molecular gas, concentrated 
in the central 175 pc (Meier & Turner 2004). The distance to NGC 6946 is estimated at 
5.5 Mpc (TuUy & Fisher 1988), to 5.9 ± 0.4 Mpc (Karachentsev, Sharina, & Huchtmeier 
2000, based on blue supergiants) ; we adopt 5.9 Mpc. 

In Figure 5 we show the radio continuum maps and F-band and Ha HST WFPC2 
images of NGC 6946. They were registered by ahgning the 6 cm peak and F-band and 
Ha peak. We identify five sources with intensity peaks in the emission complex, A - D 
and I, in the central 60 pc of NGC 6946. The most luminous 2 cm source. Source B, is 
considered to be with the strongest Ha emission. These five sources appear related to giant 
molecular cloud CMC El in Meier & Turner (2004), a cloud that has molecular gas mass of 

14 X lO^M©. Sources E, F and H are further away from the center, and are detected only at 
6 cm. They have no optical counterparts although E might be associated with the diffuse 
optical emission north-west of the optical peak. These nine compact sources contribute 
30 % of the total 6 cm emission, and 30 % of the 2 cm emission, from the central kpc of 
this galaxy. The high fraction of 6 cm emission from compact sources as compared with the 
other three spiral galaxies suggests that there are compact and young SNRs in the center 
of NGC 6946, which may have a more evolved star cluster than the other spirals. 

Sources A - D and I are within 2", or 60 pc in projection, of the peak in the HST 
F-band image. The infrared peak shown in an archived NICMOS F160W image hes on 
source I to within the pointing uncertainties of HST. Even with sub-arcsecond resolution, 
we cannot resolve the extended emission between these sources. We derive their spectral 
indices based on their peak fluxes. It is possible that they are physically connected to each 
other and that the two nonthermal sources, A and C, are the consequence of extensive star 
formation near source B, D, and I. The optical peak of NGC 6946 is associated with a 12 - 

15 Myr old Cluster 1094 (Larsen et al. 2002) of My — —10.18 mag, which is comparable in 
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luminosity to the bright open star clusters but inconspicuous in comparison with a young 
SSC. By comparison, Source B, D, and I together should have My ~ —12.5 mag. Using the 
column density of molecular hydrogen derived from CO (Meier & Turner 2004), we estimate 
that Ay ~ 20 — 30 mag toward source D and I. Source B, D, and I are probably embedded 
in CMC El. 

Sources A and C have negative spectra and are likely SNRs (Weiler et al. 1986). The 
2 cm map shows a decrease in flux density at the position of Source C. This suggests the 
contamination from extended emission might be responsible for the SNR-like spectral index 
of these two sources. Thus, we cannot distinguish them from RSNe, SNRs, or H ll regions. 
They could be explained by one or few short-lived young Type I RSN (Weiler et al. 1998). 

Source E shows extended emission at 2 cm and is marginally detected. In Table 4, we 
put its peak and integrated fluxes at 2 cm in the area of the 6 cm beam. The flat spectrum 
suggests that Source E is a large (30 Dor-like) Hll region. We do not have significant 
detections toward Sources F - H at 2 cm, and cannot identify their nature. 

6. Conclusions 

We have imaged compact radio emission in the centers of four nearby spiral galaxies, 
IC 342, Maffei II, NGC 2903, and NGC 6946, at subarcsecond resolution at 6 and 2 cm 
with the VLA. Compact radio sources contribute ~ 20 - 30% of the total 2 cm emission 
from the central kpc of these galaxies, but only 5 - 10 % at 6 cm; in NGC 6946, they 
contribute ~ 30 % at both wavelengths. The 2 cm emission within the central 100 - 200 pc 
of these spirals is dominated by H ll regions, while diffuse synchrotron emission dominates 
at 6 cm. Over half of the 38 identified compact sources appear to be Hll regions, based 
on their radio spectra; of these, one-third to half appear to have optically thick free-free 
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emission. The high resolution of these maps selects for these compact, high brightness Hll 
regions. The brightness temperatures of these optically thick, compact H ll regions are 
consistent with this identification, and indicate actual source diameters of about 1-3 pc 
(~ 0.02" - 0.05".) The compact Hll regions contain ionized gas with densities as high as 
rii ~ 10^ cm~^, which suggests that these Hll regions are relatively young. Few of the 
compact radio sources have candidate counterparts in the visual wavebands. Comparison 
with optical HST images show that the young H ll regions tend to be found in dust lanes, 
and probably are subject to high extinction and also dynamical influences from molecular 
clouds. The predicted infrared emission from these H ll regions contributes 3 - 5 % of the 
total infrared luminosity of the galaxies. The largest H ll regions require equivalent of ~ 
500 - 1,000 07 stars to excite them; since they are dense, compact Hll regions, by analogy 
with Galactic compact H ll regions, they are likely to be very youngest of the massive young 
clusters, a Myr or less in age. 
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Bigbt Ascension (J§O.Op) 

Fig. 1. — 2MASS if -band images and VLA 6 cm radio continuum maps of our four sample 
galaxies. On the left are 3.6' x 3.6' visible images from 2MASS if-band overlaid with 6 cm 
radio emission from the central region. On the right are 6 cm maps with 16" x 16" field-of- 
view, contoured at half-integral powers of 2 x 0.14 mJy beam~^ (~ 4 cr). Peak flux densities 
are 1.49, 3.13, 0.37, and 2.74 mJy/beam for IC 342, Maffei II, NGC 2903, and NGC 6946, 
respectively. Beams are ^.'3 in size. 
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Fig. 1. — Continued. 
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Fig. 2. — Closeups of the central region of IC 342. Upper-left is (a) 6 cm contours overlaid 
on a WFPC2 F-band (F555W) image. Upper-right, lower-left, and lower-right are (b) 6 cm 
and (c) 2 cm contours overlaid on a WFPC2 Ha (F656N) image. The beam size of each 
map has been matched to 0.45" x 0.42", PA= 46°. Contours are half-integral powers of ± 2 
X 0.30 and 0.43 mJy/beam for 6 cm and 2 cm, respectively. Uncertainty in the registration 
of the VLA and HST maps is ~ 1". 
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Fig. 3. — Closeups of the central region in Maffei II. Upper-left is (a) 6 cm contours overlaid 
on a WFPC2 V band (F606W) image. Upper-right and lower-left are (b) 6 cm and (c) 2 
cm contours overlaid on a WFPC2 I band (F814W) image. The beam size of each map has 
been matched to 0.43" x 0.42", PA= 40°. Contours are at half-integral powers of ± 2 x 0.21 
and 0.60 mJy/beam for 6 cm and 2 cm maps, respectively. Uncertainty in the registration 
of the VLA and HST maps is ~ 1". 
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Fig. 4. — Closeups of the central region in NGC 2903. Upper-left is (a) 6 cm contours 
overlaid on a WFPC2 F-band (F555W) image. Upper-right and lower-left are (b) 6 cm and 
(c) 2 cm contours overlaid on a WFPC2 Ha (F656N) image. The beam size of each map 
has been matched to 0.43" x 0.41", PA= 82°. Contours are half-integral powers of ± 2 x 
0.11 and 0.17 for 6 cm and 2 cm, respectively. Uncertainty in the registration of the VLA 
and HST maps is ~ 1". 
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Fig. 5. — Closeups of the center region of NGC 6946. Upper-left is (a) 6 cm contours overlaid 
on a WFPC2 F-band (F606W) image. Upper-right and the lower-left are (b) 6 cm and (c) 
2 cm contours overlaid on WFPC2 Ha (F656N) image. The orange contour represents the 
region of decreasing flux density toward the center. The beam size of each map has been 
matched to 0.42" x 0.38", PA= -79°. Contours are at half-integral powers of ± 2 x 0.10 
and 0.50 mJy/beam for 6 cm and 2 cm maps, respectively. Uncertainty in the registration 
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Fig. 6. — Luminosity function of the compact Hll regions in the centers of IC 342, Maffei 
II, NGC 6946, and NGC 2903. Filled circles with statistical error bars show the cumulated 
number counts of Hll regions found in the 4 spiral galaxies of our sample. Solid lines 
represent the best-fit luminosity function model for sources with L^a > 2 x 10'^^ erg sec^^. 
The sample is complete down to N^y^ ~ 6 x 10^° {Ha ~ 8.3 x 10^^ erg s~^). Each 07 
star contributes roug hly NLyc ~ 10^^ s-^ the brightest srouces (in Maffei II and NGC 2903) 
correspond to clusters of ~ 500 - 1,000 O stars. 
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Table 1, 


. Properties 


of the Sample 


Galaxies 






Galaxy 


Hubble Type 


(Mpc) 


Total 
(mjy) (mjy) 


'-'6 '-'2 

Mapped 
(mJy) (mJy) 


1Q'>Lq 


M{H2Y 
IO^Mq 


IC 342 

Maffci II 
NGC 2903 
NGC 6946 


SAB(rs)cd 

SAB(rs)bc 

SB(s)d 

SAB(rs)cd 


3.3 

5.0 
8.9 
5.9 


82 
107 

35 
39 


38 

46 
12 
23 


15 
16 

5 

17 


7 

13 
6 
8 


3.0 
12.3 
9.4 
7.2 


6.8 
3.7 
4.5 
4.9 



^Reference papers for distance: IC 342: Saha, Claver, & Hoessel (2002); Maffei II: Hurt, Merrill, 
Gatley, & Turner (1993); NGC 2903: Drozdovsky & Karachentsev (2000); NGC 6946: Karachentsev, 
Sharina, & Huchtmeier (2000). 

''Total radio continuum fluxes in VLA low resolution maps using B-config. at 6 cm and C-config at 
2 cm for IC 342, NGC 6946 (Turner & Ho 1983), and Maffei II (Turner & Ho 1994); using C-config. at 
6 cm for NGC 2903 (Wynn- Williams & Becklin 1985). 

'^Total fluxes in naturally- weighted, 0'.'3 maps (this paper). 

'^LjB.AS (infrared luminosity) calculated using the prescription from Sanders & Mirabel (1996). IRAS 
data are from point source flux densities listed in the IRAS Faint Source Catalog, version 2.0 (Moshir 
& et al. 1990). IRAS fluxes of Maffei II are from Chini, Kreysa, Kruegel, & Mezger (1986). 

''Molecular mass calculated using ^^CO data (Sage, Shore, & Solomon 1990) with distance correction. 



Table 2. VLA Data 



Object/Band* Obs. Date On Source Config. Flux Calib. Phase Calib. Prog. ID 

mm/dd/yyyy (sec) 

IC 342 



NGC 



. C-band 


11/23/1996 


5290 


A 


3C48 


0224-t-671 


AW419 




11/20/1988 


4040 


A 


3C48 & 3C147 


0224-1-671 


AT98'' 




07/11/1981 


6260 


A 


3C286 


0212-1-735 


AT21 


. U-band 


08/23/1999 


4310 


A 


3C286 


0228-1-673 


AT227'' 




12/08/1995 


4550 


B 


3C48 


0224+671 


AW419 


)i II 


11/14/1981 


6020 


C 


3C286 


0224-1-671 


AT21 


. C-band 


08/26/1989 


6860 


BC 


3C138 & 3C48 


0224-f671 


AH353 




11/20/1988 


4020 


A 


3C48 & 3C147 


0224-f671 


AT98'' 




09/05/1982 


2970 


B 


3C286 


0224-1-671 


AT21 




09/04/1982 


5570 


B 


3C286 


0224+671 


AT21 


. U-band 


10/12/1983 


2330 


A 


3C48 & 3C138 


0224+671 


AH14lb 




04/18/1983 


4090 


C 


3C286 


0224+671 


AT25 




04/18/1983 


1170 


c 


3C286 


0224+671 


AT25 


2903 














. C-band 


01/11/1997 


4340 


A 


3C286 


0851+202 


AW419 




11/26/1996 


4320 


A 


3C286 


0851+202 


AW419 




11/18/1988 


3470 


A 


3C286 


0851+202 


AT98^ 




02/01/1988 


870 


B 


3C286 


0953+254 


AS314 




08/16/1986 


870 


B 


3C48 


0953+254 


AD176 




04/26/1985 


5530 


B 


3C286 


0851+202 


AB324 




10/24/1982 


8550 


B 


3C286 


0851+202 


AT21 


. U-band 


08/23/1999 


4280 


A 


3C286 


0854+201 


AT227'' 




04/14/1996 


1500 


C 


3C286 & 3C147 


1004+141 


AQll 




01/05/1996 


3840 


B 


3C48 


0851+202 


AW419 




12/27/1995 


1800 


B 


3C286 


0851+202 


AW419 




12/08/1986 


17100 


C 


3C286 


0851+202 


AB324 




09/18/1985 


9480 


C 


3C286 


0851+202 


AB324 




04/18/1983 


2210 


C 


3C286 


0851+202 


AT25 




04/17/1983 


2900 


C 


3C286 


0851+202 


AT25 
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Table 2 — Continued 



Object/Band^ 


Obs. Date 
mm/dd/yyyy 


On Source 
(sec) 


Config. 


Flux Calib. 


Phase Calib. 


Prog. ID 


NGC 6946 














C-band 


06/17/1999 


3680 


AD 


3C48 


2021+614 


AS568 




11/23/1996 


4730 


A 


3C48 


2021+614 


AW419 




04/10/1994 


6340 


A 


3C286 


2021+614 


AS525 




11/20/1988 


4010 


A 


3C48 & 3C147 


2021+614 


AT98'' 


U-band 


12/08/1995 


4580 


B 


3C48 


2021+614 


AW419 




10/12/1983 


2600 


A 


3C48 & 3C138 


2021+614 


AHUl^ 




11/14/1981 


6020 


C 


3C286 


2021+614 


AT21 



'^VLA C-band: 6 cm; [/-band: 2 cm; K-hand: 1.3 cm. 
''Unpublished data from observations executed by the authors. 
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Table 3. Compact Continuum Source Observed Quantities 



Source 




a(J2000) 
C) 


(5(J2000) 
{") 


QPeak;a,c 
6 cm 

(mjy beam~^) 


C!„ b,c 

(mJy) 


QPeak;a,c 
2 cm 

(mJy beam~^) 


S2 

(mJy) 


Size<i 
// ^ // 


IC 342 




03'"46™ + 


68°05'+ 














A 


48.46 


48.4 


1.4 


1.5 ± 0.2 


0.6 


0.9 ± 0.3 


0.2 0.1, 40 




B 


47.88 


44.1 


0.8 


0.8 ± 0.2 


< 0.3« 


0.4 ± 0.3 


0.3 0.1, 120 




C 


47.73 


45.7 


0.8 


1.3 ± 0.3 


1.0 


2.2 ± 0.5 


0.4 0.3, 20 




D 


47.93 


47.3 


0.8 


1.5 ± 0.4 


< 0.3« 


< 0.3 


0.6 0.4, 10 




E 


47.61 


43.7 


0.5 


0.7 ± 0.3 


0.7 


0.7 ± 0.5 


0.4 0.3, 40 




F 


47.98 


47.8 


0.5 


0.6 ± 0.3 


< 0.3<= 


< 0.3 


0.4 0.3, 70 




G 


48.69 


44.0 


0.5 


0.6 ± 0.3 


0.5 


0.5 ± 0.5 


0.5 0.3, 100 




H 


49.18 


50.3 


0.4 


0.5 ± 0.3 


< 0.3« 


< 0.3 


0.4 0.3, 140 




I 


47.86 


45.5 


< 0.2« 


< 0.2 


0.7 


0.7 ± 0.5 


0.4 0.3, 120 




J 


48.50 


47.2 


< 0.2« 


< 0.2 


0.7 


0.7 ± 0.5 


0.8 0.2, 10 




K 


48.17 


48.5 


< 0.2= 


< 0.2 


0.7 


0.7 ± 0.5 


0.4 0.3, 150 




L 


48.99 


46.4 


< 0.2O 


< 0.2 


0.7 


0.7 ± 0.4 


0.4 0.2, 70 


Maffei II 




02''41'"+ 


59° 36'+ 














A 


55.16 


17.3 


2.6 


3.1 ± 0.2 


2.1 


2.7 ± 0.5 


0.3 0.2, 60 




B 


55.14 


21.5 


1.1 


1.1 ± 0.2 


< 0.3" 


< 0.3 


0.4 0.3, 60 




C 


55.10 


21.2 


1.0 


1.6 ± 0.2 


1.7 


2.2 ± 0.6 


0.4 0.3, 70 




D 


55.09 


18.1 


1.0 


1.5 ± 0.3 


1.3 


1.8 ± 0.7 


0.5 0.4, 150 




E 


55.17 


13.9 


0.8 


0.8 ± 0.2 


1.2 


1.3 ± 0.5 


0.3 0.2, 160 




F 


55.04 


20.5 


0.8 


1.6 ± 0.2 


1.0 


1.0 ± 0.4 


0.2 0.1, 70 




G 


55.10 


13.5 


0.8 


1.2 ± 0.2 


0.8 


1.5 ± 0.4 


0.5 0.1, 




H 


55.16 


15.5 


0.7 


1.0 ± 0.3 


0.8 


1.4 ± 0.7 


0.5 0.4, 10 




I 


55.07 


15.6 


0.6 


0.9 ± 0.3 


0.8 


0.8 ± 0.8 


0.6 0.4, 40 




J 


55.15 


15.9 


0.6 


0.7 ± 0.3 


< 0.3" 


< 0.3 


0.5 0.4, 80 


NGC 2903 




09'»32™+ 


21°30'+ 














A 


09.85 


05.2 


0.3 


0.4 ± 0.1 


0.2 


0.4 ± 0.2 


0.4 0.2, 150 




B 


10.10 


06.5 


0.3 


0.5 ± 0.2 


0.4 


1.2 ± 0.4 


0.6 0.4, 120 




C 


09.93 


-01.9 


0.3 


0.4 ± 0.2 


0.3 


0.3 ± 0.3 


0.4 0.3, 100 




D 


09.97 


06.3 


0.2 


0.2 ± 0.1 


< 0.1« 


< 0.1 


0.2 0.2, 60 
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Table 3 — Continued 



Source 




q(J2000) 

n 


<5(J2000) 
(") 


oPeak ;a,c 
6cm 

(mjy beam~^) 


(mJy) 


oPeak ;q 

(mJy beam" 


■,.c 


(mJy) 


Size'* 

// w // o 




E 


10.21 


07.9 


0.2 


0.2 ± 0.1 


0.2 




0.4 ± 0.2 


0.4 0.2, 140 




F 


09.76 


01.3 


< 0.1" 


< 0.1 


0.3 




0.4 ± 0.3 


0.5 0.4, 90 




G 


10.00 


07.5 


< 0.1« 


< 0.1 


0.2 




0.3 ± 0.2 


0.4 0.2, 50 


NGC 6946 




20*" 34'"+ 


60°09'+ 
















A 


52.34 


14.2 


1.5 




1.1 






confused 




B 


52.28 


14.2 


1.4 




1.5 






confused 




C 


52.28 


14.5 


1.1 




0.7'' 






confused 




D 


52.25 


14.7 


1.0* 




1.0 






confused 




E 


52.03 


15.0 


0.3 


0.5 ± 0.2 


0.5 




0.5 ± 0.5 


0.5 0.3, 120 




F 


52.39 


11.9 


0.2 


0.3 ± 0.1 


< 0.4<= 




< 0.4 


0.3 0.2, 70 




G 


52.15 


14.5 


0.2 


0.3 ± 0.3 


< OA" 




< 0.4 


0.9 0.5, 130 




H 


52.50 


13.6 


0.1 


0.3 ± 0.1 


< 0.4« 




< 0.4 


0.4 0.2, 140 




I 


o2,2:-! 


11.2 


0. 




1.:-! 









'*'^6cm'° and S^^^: peak flux density (intensity) of 6 cm and 2 cm continuum, respectively, with beams 0.45" X 0.42", 
PA = 46° for IC 342, 0.43" x 0.42", PA = 40° for Maffei II, 0.43" x 0.42", PA = 82° for NGC 2903, and 0.42" x 0.38", 
PA = -79° for NGC 6949. The rms noise levels at 6 cm are 0.05, 0.04, 0.03, and 0.03 mjy/beam for IC 342, Maffei 2, 
NGC 2903, and NGC 6946, respectively, and at 2 cm are 0.09, 0.11, 0.05, and 0.14 mJy/beam. 

^3(,cm and S^cm' total integrated fluxes at 6 cm and 2 cm, respectively. Quoted uncertainties in the fluxes are derived 
from the rms and measured source size; intensities are more accurate for extended sources. The systematic error of flux 
is expected to be < 5 %. 

'^Flux densities are measured in (u,ij)-restricted map. The uncertainties quoted here are 3 cr See text for details. 
•^Single Gaussian fltting uncertainty is < 30%, using AIPS task IMFIT. 
a upper limit. 

^Peak flux density is averaged flux density over a beam area. 
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Table 4. Compact Continuum Source Derived Properties 



Source 



ae-2 



Type 



Free- Free 

T 



■vrthin;6 
Lyc 



No7'= Lir" 

lO^'Lo 



Comments 



IC 342 



A 


-U.O _Q 2 


SNR 


N/A 








B 


-0 6 * 
-0.6 


SNR 


N/A 








C 


5 

U.O _(j 2 


Hli 


Thick 


2.1 


210 


4.2 


D 


< -1.5 


SNR/RSN 


N/A 








E 


0.0 1° j 


Hii ? 


Thin? 


0.7 


70 


1.4 


F 


< -0.6 


SNR 


N/A 








G 


Q o +0.6 


SNR 


N/A 








H 


< -0.5 


SNR 


N/A 








I 


> 1.1 


Hii 


Thick 


0.7 


70 


1.4 


J 


> 1.1 


Hii 


Thick 


0.7 


70 


1.4 


K 


> 1.1 


Hii 


Thick 


0.7 


70 


1.4 


L 


> 1.1 


Hii 


Thick 


0.7 


70 


1.4 


Total 








5.6 


560 


11.1 


Maflfei II . . . 














A 


-0.1 


Hii 


Thin 


6.1 


610 


12.2 


B 


< -1.2 


SNR/RSN 


N/A 








C 


U.J _Q 2 


Hii 


Thick 


4.9 


490 


9.7 


D 


2 


Hii 


Thin? 


4.1 


410 


8.2 


E 


-0.3 


Hii 


Thick? 


2.9 


290 


5.8 


F 


-0 4 


SNR 


N/A 








G 


2 

"■^ -0.2 


Hii 


Thin? 


3.4 


340 


6.8 


H 


U.CS _Q 4 


Hii 


Thick? 


3.2 


320 


6.5 


I 


-U.i _o.7 


Hii ? 


Thin? 


1.9 


190 


3.7 


J 


< -0.7 


SNR 


N/A 








Total 








26 


2600 


53 


NGC 2903 . . . 














A 


0.0 t°oi 


Hii ? 


Thin? 


2.9 


290 


5.8 


B 


0.8 1° j 


Hii 


Thick 


8.2 


820 


16 


C 


-0 2 +°-* 


Hii 


Thin? 


2.3 


230 


4.6 


D 


< -0.5 


SNR 


N/A 
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Table 4 — Continued 



Source 


0:6-2'' 


Type 


Free- Free 
r 


Tvrthin;6 

lOSi s-1 


No7'= 


LiRd 

lO'-Lo 


Comments 


E 
F 


0.6 1^4 
> 1.2 


Hii 
Hii 


Thick 
Thick 


2.9 
3.1 


290 
310 


5.8 
6.1 




Q 


^ 1 n 


H II 


Thick 


2.1 


210 


4.2 




Total 








21 


2100 


43 




NGC 6946 . . . 
















A« 

B" 
C 


-0 3 +°-^ 

-U.J _g 2 

1 

"■^ -0.2 
< -0.4 


SNR? 
Hii ? 
SNR? 


N/A 
Thin? 
N/A 


4.5 


450 


9.0 


hole at 2cm 


D« 

E 

F 




U.U _Q 2 

+° '' 
-0.4 

< 0.3 


Hii ? 
Hii ? 
? 


Thin? 
Thin? 


3.0 
1.5 


300 
150 


6.0 
2.7 


undetermined 


G 


< 0.3 


? 










undetermined 


H 


< 0.3 


? 










undetermined 


r 


> 1.1 


Hii 


Thick 


4.5 


450 


9.1 




Total 








14 


1400 


27 





"Spectral indices are computed from the total flux densities. The upper and lower limits are based 
on \/2 a release at both wavelengths. 

''Assume the electron temperature Te = 10,000 K. The values will be ~ 15% lower than listed if 
Te = 6,000 K is assumed. 

"^Assume the ionization rate of a reference O star. Nor ^ Nz^yc/W^^ sec~^ . 

Using Ljpi/NLye = 2x 10~''''(Lq sec) for Lju estimate, based on ZAMS cluster (following Turner, 
Ho, & Beck 1998), and assuming Lm = Lqb- 

''confused sources, a are derived from peak fluxes. The Nj^", Nor, a^id ^/H of them might be 
overestimated due to possible contamination in peaJc flux measurements by extended emission. 



